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Abstract 

A  Static  extended  trailing  edge  attached  to  a  NACA0012  airfoil  section  is  studied  for  achieving  lift  enhancement  at  a 
small  drag  penalty.  It  is  indicated  that  the  thin  extended  trailing  edge  can  enhance  the  lift  while  the  zero-lift  drag  is 
not  significantly  increased.  Experiments  and  calculations  are  conducted  to  compare  the  aerodynamic  characteristics 
of  the  extended  trailing  edge  with  those  of  Gurney  flap  and  conventional  flap.  The  extended  trailing  edge,  as  a 
simple  mechanical  device  added  on  a  wing  without  altering  the  basic  configuration,  has  a  good  potential  to  improve 
the  cruise  flight  efficiency. 
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Introduction 

It  has  long  been  recognized  that  flaps  can  significantly  alter  wing  aerodynamics  for  high  lift  generation. 
Conventional  flaps  have  larger  drag;  therefore  they  are  mainly  deployed  for  takeoff  and  landing  and  are  not  suitable 
in  cruise  flight.  Small  trailing  edge  devices  like  Gurney  flaps  and  divergent  trailing  edges  for  lift  enhancement  have 
attracted  considerable  attention  in  aeronautical  community.  The  Gurney  flap  was  introduced  by  Liebeck  [1]  for 
aeronautical  applications,  and  considerable  measurements  and  calculations  have  been  performed  to  determine  the 
aerodynamic  characteristics  of  wings  with  Gurney  flaps  at  low  speeds  [2-9].  Measurements  were  also  made  at  high 
speeds  [10].  Gurney-flap-type  tabs  have  been  used  along  with  vortex  generators  to  enhance  lift  and  suppress  flow 
separation  for  a  single- slotted-flap  airfoil  at  large  flap  deflection  [11].  Microflaps  that  can  change  the  deflection 
angle  were  used  for  flutter  suppression  [12].  To  mimic  the  behavior  of  bird  feathers  during  landing,  a  self-adjusting 
movable  flap  close  to  the  trailing  edge  was  proposed  by  Bechert  et  al.  [13],  which  was  able  to  maintain  higher  lift 
when  stall  occurs.  A  review  on  trailing  edge  devices  and  other  control  technologies  is  given  by  Stanewsky  [14]. 

Measurements  of  avian  wing  geometry  for  extracting  typical  avian  airfoil  sections  indicate  that  the  merganser  and 
owl  wings  have  a  very  thin  trailing  edge  that  is  a  single  layer  of  feathers  extended  from  a  ‘normal’  airfoil  section 
[15].  Naturally,  a  question  is  whether  and  how  this  unique  trailing  edge  affects  the  global  aerodynamic 
characteristics  of  the  wing.  This  inspires  the  concept  of  a  static  extended  trailing  edge  (SETE)  as  illustrated  in 
Eigure  1(a).  The  main  airfoil  element  remains  unchanged,  but  it  is  extended  at  the  trailing  edge  by  attaching  a  thin 
splitting  plate  of  suitable  length  and  rigidity.  As  shown  in  Eig.  1(b),  the  basic  geometrical  parameters  are  the  main 
airfoil  chord  (c),  length  of  SETE  (/)  and  deflection  angle  ( ^^ ).  Depending  on  specific  applications,  the  extension 
could  be  an  aluminum  plate,  polymer  membrane,  composite  sheet  or  smart  material  plate.  In  general,  a  thin 
extended  trailing  edge  is  flexible  such  that  it  can  be  passively  changed  through  flow- structure  interaction  or  actively 
controlled  by  embedded  actuators.  In  the  current  operating  conditions,  our  analysis  shows  that  when  extended 
trailing  edges  are  sufficiently  rigid,  their  shapes  are  not  sufficiently  affected  by  flow  and  no  vibration  is  induced. 
Thus,  the  cases  studied  here  are  static  or  quasi-steady.  The  effects  of  SETE  on  the  wing  aerodynamics  are  mainly 
due  to  modifications  of  the  airfoil  camber  and  of  the  flow  structure  at  the  trailing  edge.  The  lift  enhancement  by 
SETE  is  expectedly  due  to  the  camber  effect.  It  is  speculated  that  the  drag  penalty  of  SETE  is  small  since  it  is 
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usually  embedded  in  the  wake  of  the  main  airfoil.  The  drag  penalty  should  be  examined  carefully  in  comparison 
with  Gurney  flaps  and  conventional  flaps.  The  objective  of  this  work  is  to  study  SETE  as  an  alternative  flow  control 
technique  for  lift  enhancement  at  a  small  drag  penalty  in  cruise  flight. 

We  first  describe  the  geometrical  features  of  SETE  and  relevant  similarity  parameters.  Experimental  results  are 
presented  including  the  lift  and  drag  coefficients  of  the  baseline  NACA0012  airfoil  model  and  model  with  SETE  and 
Gurney  flaps.  The  lift  enhancement  is  achieved  by  SETE,  while  the  drag  polar  and  lift-to-drag  ratio  remains  largely 
unchanged.  The  zero-lift  drag  and  Oswald  efficiency  are  not  significantly  affected  by  SETE.  The  benefit  analysis 
for  SETE  in  cruise  flight  is  given  in  comparison  with  Gurney  flap.  Based  on  thin  airfoil  theory  and  CED 
simulations,  the  aerodynamic  characteristics  of  SETE  are  further  compared  with  those  of  the  corresponding  Gurney 
flap  and  conventional  flap. 


(a)  (b) 

Figure  1.  (a)  Illustration  of  an  extended  trailing  edge,  (b)  Geometrical  parameters  and  coordinate  system. 

The  team  members  of  experimentalists  and  computationalists  are  fully  engaged  throughout  the  various  phases  of  the 
research.  The  open  and  continuous  communications  have  facilitated  the  planning,  the  design,  and  the  execution  of 
experimental  and  computational  tasks.  In  fact,  they  often  lead  to  valuable  insights  that  are  not  readily  apparent  to  a 
uni-approach  team.  We  believe  a  full  integration  of  experimental  and  numerical  approaches  can  be  very  effective 
and  cost-efficient  for  aerodynamics  studies. 

Wind  Tunnel  and  Model 

Tests  were  conducted  in  the  Advanced  Design  Wind  Tunnel  (ADWT)  at  Western  Michigan  University  (WMU). 
The  ADWT  is  a  low  speed,  closed  circuit  single  return,  atmospheric  tunnel  with  a  test  section  of  32  inches  high,  45 
inches  wide,  and  8  feet  long.  Eull  test  section  length  Plexiglas  windows  on  sides  and  ceiling  allow  optical  access. 
The  161  feet  long  circuit  consists  of  two  legs  of  65.6  feet  long  and  the  other  two  legs  of  14.9  feet  long.  Low 
turbulence  and  good  flow  direction  control  are  obtained  by  a  combination  of  features,  including  four  screens  and 
honeycomb  in  the  stilling  camber,  a  10:1  contraction  ratio  between  the  stilling  chamber  and  test  section,  empirically 
tailored  test  section  walls,  tangential  flow  breathing  at  the  test  section  exit  to  the  diffuser,  and  a  5.4°  diffuser  angle. 

Under  the  present  test  conditions  at  Re^  -4.74xl0\  natural  boundary  layer  transition  was  observed  at  about  30%c 
on  the  NACA0012  airfoil  model  at  zero  AoA  based  on  global  luminescent  oil  film  skin  friction  measurements  [16]. 
The  main  drive  power  source  is  a  125  hp  DC  motor  with  a  solid-state  speed  controller.  The  16-blade  fan  is  6  feet  in 
diameter.  The  flow  velocity  range  is  approximately  20  to  240  ft/sec  (6  to  73  m/sec)  where  the  total  temperature  is  in 
the  range  of  60  to  120°  E.  This  corresponds  to  the  dynamic  pressures  up  to  65  Ib/ft^  and  Reynolds  numbers  up  to  1.3 
X  10^  per  foot. 

A  six-component  internal  force  balance  (AeroLab)  is  installed  in  the  ADWT.  The  balance  measurement  ranges  are 
±150  lb  for  lift,  ±100  lb  for  side  force,  ±25  lb  for  axial  force,  ±250  in-lb  for  pitching  moment,  ±180  in-lb  for  yawing 
moment,  and  ±200  in-lb  for  rolling  moment.  A  parallelogram  linkage  sting  support  system  is  able  to  set  the  angle  of 
attack  from  -36°  to  36°  within  the  accuracy  of  0.1°.  A  traverse  system  of  two  degrees  of  freedom  is  equipped  for 
probes. 

A  ceramic  NACA0012  airfoil  section  model  was  used  in  experiments.  The  chord  and  span  of  the  model  were  10  in 
and  12  in,  respectively.  Along  the  wing  span,  two  1.875-in  long  and  6-in  wide  parts  made  of  polycarbonate  can  be 
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mounted  to  the  baseline  model  with  4  screws  for  clamping  SETE  at  the  trailing  edge,  and  they  can  be  easily 
removed  to  place  different  extended  trailing  edges.  Several  pieces  of  SETE  with  different  lengths  were  prepared  for 
tests,  which  were  made  of  either  aluminum  or  Mylar  (PET-polyester).  The  Young’s  modulus  and  Poisson  ratio  for 
aluminum  were  70  GPa  and  0.35,  respectively,  while  the  Young’s  modulus  and  Poisson  ratio  for  Mylar  were  2.8 
GPa  and  0.37,  respectively.  Eigure  2  shows  the  NACA0012  airfoil  model  attached  with  a  Mylar  SETE,  and  the 
model  parts  before  assembly,  including  the  baseline  NACA0012  model,  two  removable  parts  for  clamping  SETE, 
and  typical  aluminum  SETE  and  Mylar  SETE.  In  order  to  support  the  model,  the  end  plates  (18  by  13  in)  were 
mounted  to  the  sting  balance  through  two  aluminum  bars  placed  in  a  mounting  bracket.  The  end  plates  also  reduced 
the  three-dimensionality  of  flow.  Eigure  3  shows  the  model  mounted  on  the  balance  for  testing  in  the  ADWT. 


(a)  (b) 

Figure  2.  (a)  NACA0012  airfoil  model  with  a  Mylar  SETE,  (b)  Baseline  model,  removable  parts,  and 
aluminum  SETE  and  Mylar  SETE. 


Figure  3.  NACA0012  airfoil  model  mounted  on  the  end  plates  connected  to  in  a  sting  balance. 

Similarity  Parameters 

The  similarity  parameters  for  a  rigid  airfoil  with  an  extended  trailing  edge  are  considered.  The  similarity  parameters 
for  a  rigid  airfoil  are  Reynolds  number  Re^  -  U^c  /  v  and  Mach  number  ,  where  is  the  freestream  velocity 
and  c  is  the  airfoil  chord.  Since  an  extended  trailing  edge  is  attached  to  the  rigid  airfoil,  additional  similarity 
parameters  should  be  obtained  for  achieving  the  complete  flow  and  structure  similarity.  An  extended  trailing  edge  is 
generally  considered  to  be  flexible  as  a  thin  cantilever  plate  deforming  under  aerodynamic  loading.  The  differential 
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equation  for  the  time-independent  displacement  w  of  a  thin  plate  is  generally  expressed  by  L[w(  x,y  )]  =  F(  x,y ) 
[17],  where  x  and  y  are  the  chord  wise  and  normal  coordinates,  respectively,  L  is  a  linear  differential  operator,  and 
F(x,y)  is  an  external  distributed  force  on  the  plate.  In  our  case,  F( x,y )-  Ap  is  the  aerodynamic  pressure 

difference  across  the  plate  or  aerodynamic  load.  For  a  homogenous  plate,  the  operator  L  =  is  the  bi-harmonic 

operator,  where  =  Eh^  / 12(1- )  is  the  plate  flexural  rigidity.  The  rigidity  depends  on  the  plate  thickness 
(h).  Young’s  modulus  (£'),  and  Poisson  ratio  ( v). 

Let  us  introduce  the  following  the  non-dimensional  variables  F  -  F / w  =  w/ 1  ,  x  =  x/ 1 ,  y  =  y /I ,  where 
^00  ^  PMI  is  the  dynamical  pressure  and  I  is  the  length  of  SETE.  Thus,  we  have  the  non-dimensional  equation 
GjV^w  =  F  ,  where  the  main  similarity  parameter  for  aeroelasticity  is  Gj  =  /( ) .  The  functional  relations  for 

the  lift  and  drag  coefficients  are  C^  =  fj(  a,  Re^,M^,G^,  5,1/  c)  and  Cj^=  f2(  a,  Gj,S,l/c ).  For  SETE, 

the  similarity  parameter  Gj  should  be  sufficiently  large  such  that  the  shape  of  SETE  is  not  changed  significantly  by 
flow.  For  the  present  low-speed  experiments,  is  very  small  and  Re^  is  fixed  at  4.74x10^ .  Although  the 
effects  of  Re^  and  are  worthwhile  for  further  investigation,  we  mainly  consider  here  the  aerodynamic  behavior 
of  the  NACA0012  airfoil  model  with  SETE  in  the  parametric  space  (  a,  5,1/  c)  . 

Results  and  Discussions 

Static  Extended  Trailing  Edge 

Figure  4  shows  the  lift  and  drag  coefficients  as  a  function  of  angle  of  attack  (AoA)  for  the  baseline  NACA0012 
airfoil  model  at  Re^  =4.74x10^  According  to  McCormick’s  formula  [18],  the  lift  slope  depends  on  the  aspect 
ratio  (AR)  by 

a=l^  = _ ^2^ _  (1) 

da  AR  +  2(AR  +  4)/(AR  +  2)’ 

where  Gq  =  2 tt  according  to  thin  airfoil  theory.  As  shown  in  Fig.  4,  a  fit  to  the  lift  data  of  the  baseline  model  using 
Eq.  (1)  gives  the  effective  aspect  ratio  AR  =  4.42  that  is  larger  than  the  physical  aspect  ratio  of  1.2.  The  end  plates 
substantially  reduced  the  three-dimensionality  of  flow.  As  shown  in  Fig.  5,  the  drag  coefficient  is  represented  as  a 
function  of  the  lift  coefficient  squared  and  fitted  using  the  classical  relation  ^C\/  ne  AR  .  Therefore,  we 

obtain  the  zero-lift  drag  coefficient  =  0.014  and  the  Oswald  span  efficiency  e  =  0.8  for  AR  =  4.42  .  Although 
the  finite  end  plates  do  not  completely  eliminate  the  three-dimensionality  of  flow,  they  modify  the  tip  vortices  acting 
like  large  winglets. 

In  tests,  an  aluminum  sheet  (0.216  mm  thick)  and  a  Mylar  sheet  (0.254  mm  thick)  have  been  used  as  SETE,  whose 
deflection  angle  was  set  a  priori.  The  plate  flexural  rigidities  of  the  aluminum  sheet  and  Mylar  sheet  are 
=0.067  Nm  and  =0.0044  Nm  ,  respectively.  Figure  6  shows  the  lift  coefficient  as  a  function  of  AoA  and 

drag  polar  for  the  model  with  SETE  at  Re^  =4.74xl0\  where  Q  and  are  based  on  the  projected  wing  area  of 
the  baseline  model.  The  lengths  of  SETE  are  5%  and  10%  of  the  chord  of  the  baseline  model,  and  the  preset 
deflection  angles  of  SETE  are  0,  5,  7  and  10  degrees.  According  to  Fig.  6,  the  Q  distribution  is  shifted  up  and  the 
lift  is  enhanced  depending  on  the  deflection  angle  and  relative  length  of  SETE,  while  the  drag  polar  basically 
remains  unchanged.  The  drag  polar  curves  are  collapsed  for  different  deflection  angles,  indicating  that  the  zero-lift 
drag  and  Oswald  efficiency  are  not  changed.  Stall  occurs  a  few  degrees  earlier  particularly  for  larger  deflections. 
For  the  model  with  zero  deflection,  Q  is  almost  the  same  as  that  of  the  baseline  model  when  AoA  is  less  than  10°, 
but  it  is  larger  for  higher  AoA  until  stall.  Typically,  four  and  five  repeated  tests  were  conducted  for  each  case,  and 
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the  results  were  averaged.  The  measurement  uncertainties  in  Q  are  0.001  for  low  AoA  and  0.01  for  high  AoA. 
For  ,  the  measurement  errors  are  0.0002  for  low  AoA  and  0.003  for  high  AoA. 

Further,  the  drag  coefficient  is  plotted  as  a  function  of  the  lift  coefficient  squared,  as  shown  in  Fig.  7(a).  All  the  data 
are  basically  collapsed  into  a  single  straight  line,  which  is  consistent  with  Fig.  6(b).  Both  the  zero-lift  drag 
coefficient  and  Oswald  efficiency  do  not  significantly  change  within  the  measurement  accuracy  for  different 
deflection  angles  (small  enough)  and  lengths  of  SETE.  It  is  indicated  in  Fig.  7(a)  that  the  increased  results 

mainly  from  an  increase  of  the  induced  drag  C\/ neAR  associated  with  the  lift  enhancement  by  SETE.  It  is 
expected  that  this  drag  penalty  will  be  reduced  as  the  aspect  ratio  increases.  To  examine  the  aerodynamic 
effectiveness  of  SETE,  the  lift-to-drag  ratio  ( L/Z) )  is  plotted  in  Fig.  7(b)  as  a  function  of  Q  .  The  plots  of  L/ D 

against  Q  are  approximately  collapsed.  One  implication  is  that  SETE  is  able  to  increase  Q  while  the  variation  in 
L/Z)  is  relatively  small.  The  significance  of  these  findings  is  that  SETE  could  improve  the  aerodynamic  efficiency 
of  cruise  flight. 
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Figure  4.  Lift  and  drag  coefficients  as  a  function  of 
AoA  for  the  baseline  NACA0012  airfoil  model  at 

Re^  =4.74xl0\ 
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Figure  5.  Drag  coefficient  as  a  function  of  the  lift 
coefficient  squared  for  the  baseline  NACA0012  airfoil 
model  at  Re,  -4.74xl0\ 
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(a)  (b) 

Figure  6.  (a)  Lift  coefficient  as  a  function  of  AoA  and  (b)  drag  polar  for  the  NACA0012  airfoil  model  with 
SETEat  Re^=4J4xlO\ 


(a)  (b) 

Figure  7.  (a)  Drag  coefficient  as  a  function  of  the  lift  coefficient  squared,  (b)  Lift-to-drag  ratio  and  as  a 
function  of  the  lift  coefficient  for  the  NACA0012  airfoil  model  with  SETE  at  Re^  =4.74x10^ . 

Comparison  with  Gurney  Flaps 

The  aerodynamic  characteristics  of  SETE  are  compared  with  those  of  Gurney  flap,  which  is  a  short  wall  attached  to 
the  trailing  edge  that  is  perpendicular  to  the  chord  line  on  the  pressure  side  of  an  airfoil,  as  illustrated  in  Eig.  8.  Both 
SETE  and  Gurney  flap  increase  the  lift  at  the  cost  of  increasing  the  drag.  However,  there  are  some  important 
quantitative  differences  between  SETE  and  Gurney  flaps  on  the  drag  penalty.  The  measured  results  by  Li  et  al.  17] 
for  a  NACA0012  rectangular  wing  indicate  that  L/ D  as  a  function  of  Q  for  the  wing  with  Gurney  flaps  (0.5-3%c 

heights)  are  significantly  lower  than  that  for  the  baseline  wing  except  at  higher  Q .  Unlike  SETE,  the  curves  of 
L/ D  as  a  function  of  Q  for  Gurney  flaps  do  not  collapsed.  Similar  results  can  be  found  for  other  wings  with 
Gurney  flaps.  The  general  conclusion  is  that  Gurney  flaps  are  not  suitable  for  cruise  flight  due  to  the  reduced  L/ D 
by  the  larger  drag  penalty.  After  fitting  the  data  of  Li  et  al.  17]  using  +  Cf  / Tie  AR  ,  it  is  found  that  the 

zero-lift  drag  coefficient  increases  in  a  roughly  linear  fashion  from  0.0065  for  the  baseline  model  to  0.021  for 
the  3%c  high  Gurney  flap,  where  the  effective  aspect  ratio  is  46.5.  The  drag  increase  is  caused  by  the  extra  form 
drag  associated  with  open  flow  separation  behind  a  Gurney  flap. 

To  provide  a  direct  comparison  of  SETE  with  Gurney  flaps,  Gurney  flap  measurements  were  conducted  using  the 
same  baseline  NACA0012  airfoil  model  with  the  same  chord  Reynolds  number.  The  Gurney  flap  height  ranges 
from  1.2%c  to  6.75%c.  Eigure  8(b)  shows  the  NACA0012  airfoil  model  with  a  Gurney  flap.  Eigure  9  shows  the  lift 

coefficient  as  a  function  of  AoA  and  drag  polar  for  the  model  with  Gurney  flaps  at  Re^  =4.74x10^  The  Q  curve 
is  shifted  up  depending  on  the  Gurney  flap  height,  but  increases  considerably  except  for  the  cases  with  1.2%c 
and  2.0%c.  The  drag  coefficient  is  plotted  as  a  function  of  the  lift  coefficient  squared,  as  shown  in  Eig.  10(a).  The 
plots  for  the  heights  1.2%c  and  2%c  are  close  to  that  for  the  baseline  model,  while  other  curves  are  considerably 
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shifted  and  the  slope  is  also  changed.  This  indicates  that  the  zero-lift  drag  coefficient  increases  and  Oswald 
efficiency  varies  for  Gurney  flaps.  These  results  are  similar  to  the  previously  published  data  for  a  NACA0012 
airfoil  [7].  The  lift-to-drag  ratio  is  plotted  in  Fig.  10(b)  as  a  function  of  Q  .  Unlike  SETE,  the  L/ D  curves  are  not 
collapsed,  and  they  are  lower  than  the  baseline  data.  Figure  11  shows  a  direct  comparison  between  a  10%c  SETE 
with  a  5°  deflection  angle  and  a  1.2%c  Gurney  flap.  The  vertical  displacement  of  the  tip  of  the  SETE  is  about 
0.87%c  that  is  close  to  the  Gurney  flap  height  1.2%c.  The  lift  enhancements  in  the  both  cases  are  almost  the  same, 
but  the  Gurney  flap  suffers  from  a  larger  drag  penalty,  giving  a  smaller  L/ D  than  that  for  the  SETE  for  all  positive 

Q- 


To  evaluate  the  net  benefit  of  flow  control  in  terms  of  the  power  required  for  cruise  flight,  Liu  [19]  has  proposed  the 
following  weight  criterion 
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where  AAW  /W  is  the  ratio  between  the  control  system  weight  ZllT  and  the  total  aircraft  weight  W  ,  and 
/AW  is  the  actuating  power  density.  For  an  ideal,  weightless  passive  control  device,  the  criterion  is 


reduced  to 
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The  function  g  is  interpreted  as  a  benefit  margin  for  flow  control.  When  g  is  positive,  less  engine  power  is 


required  for  cruise  flight  with  a  passive  flow  control  device  like  SETE  and  Gurney  flap,  and  therefore  the  net  benefit 
is  achieved.  Figure  12  shows  the  benefit  margin  for  SETE  and  Gurney  flaps.  For  the  SETE  of  the  10%c  length  and 
deflection  angle  of  5®,  its  benefit  margin  remains  positive  and  approaches  to  zero  as  AoA  reaches  to  13°.  This 
further  confirms  the  feasibility  of  this  SETE  for  improving  the  cruise  flight  efficiency.  In  contrast,  for  the  10%c 
SETE  with  a  deflection  of  14°,  the  benefit  margin  becomes  negative  at  about  5°  of  AoA.  The  Gurney  flap  with  the 
1.2%c  height  has  a  similar  behavior,  while  for  the  Gurney  flap  with  the  4.12%c  height  g  is  considerably  below  zero 
after  about  3°  of  AoA. 
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(a)  (b) 

Figure  9.  (a)  Lift  coefficient  and  (b)  drag  coefficient  as  a  function  of  AoA  for  the  NACA0012  airfoil  model 
with  Gurney  flaps  at  Re^  =  4.74  . 
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(a)  (b) 

Figure  10.  (a)  Drag  coefficient  as  a  function  of  the  lift  coefficient  squared,  (b)  Lift-to-drag  ratio  as  a  function 
of  the  lift  coefficient  for  the  NACA0012  airfoil  model  with  Gurney  flaps  at  Re^  -4.74xl0\ 
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(a)  (b) 

Figure  11.  Comparison  between  SETE  and  Gurney  flap  on  (a)  the  lift  coefficient  as  a  function  of  AoA,  and  (b) 
the  lift-to-drag  ratio  as  a  function  of  the  lift  coefficient  squared  for  the  NACA0012  airfoil  model  at 

Re^  =4.74xl0\ 


Figure  12.  Benefit  margins  of  SETE  and  Gurney  flaps  for  cruise  flight. 


Thin  Airfoil  Theory 


Like  conventional  flaps,  the  lift  enhancement  by  SETE  is  generated  by  a  change  of  the  camber  at  the  trailing  edge. 
Eor  geometrical  simplicity,  the  main  wing  and  SETE  are  modeled  by  straight  segments,  as  shown  in  Eig.  1(b).  The 
total  projected  chord  is  c'  =  c-\-lcosS  .  The  slopes  of  the  chord  line  segments  are  0  for  0<  x<c  and  -tan  5  for 
c<x<c' .  According  to  the  classical  thin  airfoil  theory  assuming  that  the  slope  is  small  [20],  the  strength  of  the 
vortex  sheet  is 


r(0)=2u^ 


,  1  +  cos  6 

\  . .  + 

sin  6 


Sa. 


sin  nO 


(4) 


where  the  angular  variable  is  related  the  main  chordwise  coordinate  by  x  =  (c'  /  2  )(1-  cos  0 ) .  The  coefficients  in 


Eq.  (4)  are  Ag=  a  +  tanS(  1-0^/ 7i )  and  \=(2 / 7in)tand sin(n0^  )  {n  =  1,2,-').  Here  the  angular  variable  6^, 
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which  corresponds  to  the  breaking  point  between  the  main  wing  and  SETE,  is  given  by  6^  =  cos  s  - 1 ) /( s  + 1 

where  s  =  (l/c )cos 5  . 


The  sectional  lift  coefficient  is  given  by 

^ ^  -2n{l^  £)\a^  tan  S{l-0^/  7i  +  sin  0^/  7i)\.  (5) 

(l/2)pJJ^c 

The  positive  deflection  angle  of  SETE  shifts  the  Q  curve  upward  (lift  enhancement)  and  slightly  increases  the  lift 
slope.  The  moment  coefficient  around  the  1/4  total  chord  is 

Cm c'/4  = - ^ ,  =(l/4)(l  +  sftanS  (sin 26^-2 sin 6^ ) ,  (6) 

(l/2)pjjy 

where  defined  to  be  positive  when  it  makes  the  leading  edge  up.  Clearly,  the  positive  deflection  angle 

produces  a  negative  moment.  Eurthermore,  the  pressure  coefficient  difference  between  the  pressure  and  upper 
surfaces  of  the  airfoil  is  AC  ^  =2y(  x  )/U^  . 


Eor  a  finite  rectangular  wing,  the  effective  angle  of  attack  is  a-(l  +  T  )C^  / n  AR  ,  where  AR  is  the  effective  aspect 
ratio  and  the  parameter  r  is  a  function  of  AR/ 2 ti  for  a  rectangular  wing  [21].  A  convenient  regression  formula  is 
T  =  -0.0476(  AR/ 271  f  +0.2195(  AR/ 271 )  for  AR/ 271  g  [ 0.25,1.75 ]  .  Hence,  the  lift  enhancement 

AC^  =C^-  for  a  finite  rectangular  wing  can  be  calculated,  where  the  lift  coefficient  for  the  wing  with 

SETE  and  baseline  wing  are,  respectively. 
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According  to  the  above  relations,  the  lift  enhancement  by  SETE  is 
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The  difference  of  the  pressure  coefficient  between  the  pressure  and  upper  surfaces  of  the  main  wing  is  sensitively 
affected  by  SETE.  These  results  indicate  a  direct  correlation  of  the  aerodynamic  loading  between  the  main  airfoil 
and  SETE.  Eigure  13  is  a  comparison  between  the  predicted  AC^  by  thin  airfoil  theory  and  experimental  and  CED 

data  for  the  NACA0012  airfoil  section  ( AR  =  4.42  ,  l/c  =  0.1)  with  SETE.  The  near-linear  relation  given  by  thin 
airfoil  theory  reasonably  describes  the  initial  trend  of  the  experimental  results  for  small  deflections.  However,  the 
experimental  data  saturate  as  the  deflection  angle  increases  further,  which  may  be  caused  by  a  reduced  effective 
deflection  due  to  viscous  effect.  Two-dimensional  CED  calculations  using  a  Navier- Stokes  code  (Eluent),  assuming 
that  the  boundary  layer  on  the  NACA0012  airfoil  with  SETE  is  fully  turbulent,  predict  the  smaller  lift  enhancement. 


Thin  airfoil  theory  has  been  expanded  for  the  lift  problem  of  an  airfoil  with  a  Gurney  flap  [22].  A  solution  of  the 
generalized  thin-airfoil  integral  equation  for  the  vortex  strength  can  be  obtained  by  the  method  of  successive 
approximations.  The  Gurney  flap  lift  enhancement  is  interpreted  as  a  special  camber  effect.  The  lift  coefficient 
increment  (pitching  moment  as  well)  is  given  as  a  square-root  function  of  the  relative  Gurney  flap  height,  i.e., 

AC^  oc  -Jh/c  .  Eor  comparison  with  Gurney  flap,  the  vertical  displacement  of  SETE  is  defined  as 
h/c  =  (l/c)sinS.  Therefore,  AC^  by  SETE  depends  on  the  relative  length  / / c  and  deflection  angle  5  .  Eigure 
14  shows  the  normalized  AC^  as  a  function  of  h/ c  for  SETE  and  Gurney  flap  according  to  thin  airfoil  theory. 
Interestingly,  AC^  increases  in  a  linear  fashion  as  an  increase  of  the  deflection  angle  for  SETE  with  a  fixed  length. 
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while  ZlQ  follows  a  square-root  relation  for  both  Gurney  flap  and  SETE  with  a  variable  length  for  a  fixed 
deflection  angle. 


Figure  13.  Measured  lift  enhancement  compared  with 
thin  airfoil  theory  and  CFD  data  for  the  NACA0012 
airfoil  model  with  SETE. 
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Figure  14.  Comparison  of  normalized  lift 
enhancement  between  SETE  and  Gurney  flap  based 
on  thin  airfoil  theory. 


Deformation  and  Bending  Moment 

The  aluminum  and  Mylar  sheets  as  SETE  are  usually  set  at  nominal  deflection  angles  in  experiments  assuming  they 
are  sufficiently  rigid.  Nevertheless,  small  aeroelastic  deformation  relative  to  the  neutral  position  always  exists.  To 
estimate  the  elastic  deformation  of  SETE,  we  consider  a  2D  clamped  thin  plate  deformed  due  to  aerodynamic  load. 
The  non-dimensional  deformation  w  -w/ 1  is  given  by 

_  _  r  ^3  r^s  _  _  _ 

w(  X  )  =  Gj  \  dxji  dx2\  dxA  ACp(  ^4  )dx^ Bx  / 2  ,  (10) 

Jo  Jo  Jo  Jo 

where  x-x/l  is  the  normalized  chordwise  coordinate  from  the  breaking  point  between  the  main  element  and 
SETE,  G,  =  /( ) ,  A  =  2G;’ (  g,  -  g,  )  ,  A  =  Gj’ (  g,-2g,  )  ,  and 

pi  p^7  p^2  _ 

=  dxj\  dx2  AC  (  x^  )dx^  , 

Jo  Jo  Jo 

g2=\  cE^\  AC  J  X2  )dX2  . 

Jo  Jo 

The  pressure  coefficient  difference  =  2/(  x )/U^  from  thin  airfoil  theory  is  substituted  into  Eq.  (10)  and 
numerical  integration  is  carried  out.  Eigure  15  shows  the  normalized  deformation  w/  sin  S  as  a  function  of 
X  =  x/l  for  the  0.216  mm  thick  aluminum  sheet  (Gj  =5.74  ,  10%c)  and  0.254  mm  thick  Mylar  sheet  {Gj=  0.38  , 
10%c).  The  deformation  of  the  aluminum  SETE  is  less  than  1%  of  the  vertical  tip  displacement  IsinS ,  while  for 
the  Mylar  SETE  it  could  be  13%.  To  estimate  the  actuating  power  required  for  deflecting  SETE,  the  bending 
moment  relative  to  the  breaking  point  between  the  main  element  and  SETE  is  estimated.  The  bending  moment 
coefficient  is  defined  as 

M  _  _ 

Cm,  =^=  AC^(x)xdx. 
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Figure  16  shows  the  bending  moment  coefficient  as  a  function  of  AoA  for  SETE  with  a  deflection  angle  of  5 
degrees.  The  bending  moment  for  the  2D  case  is  in  the  order  of  qj^  . 


Figure  15.  Normalized  deformation  of  the  aluminum  SETE  and  Mylar  SETE  with  a  deflection  angle  of  5 
degrees  for  different  AoAs. 


Figure  16.  Bending  moment  of  SETE  with  a  deflection  angle  of  5  degrees  relative  the  breaking  point  between 
the  main  element  and  SETE. 

CFD  Simulations 

In  order  to  provide  further  insights  to  the  lift  enhancement  by  SETE,  the  flow  characteristics  around  a  NACA0012 
airfoil  with  SETE  were  studied  using  CFD  simulations  Re^  -4.74x10^ .  The  problem  of  a  steady  turbulent  flow 
over  the  airfoil  was  solved  using  the  connnercial  code  Fluent  on  a  C-grid.  Grid  independence  was  validated  by 
comparing  pressure  distributions  along  the  airfoil  at  two  different  angles  of  attack  with  published  results  [23].  The 
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lift  and  drag  coefficients  were  also  compared  with  the  published  data  [23].  The  boundary  layers  were  assumed  to  be 
fully  turbulent,  and  the  realizable  k-epsilon  turbulence  model  with  the  non-equilibrium  wall  function  was  used. 

Figure  17  shows  how  the  defection  angle  of  SETE  affects  the  pressure  coefficient  distribution  on  the  airfoil  for  zero 
AoA.  The  aerodynamic  loading  on  SETE  itself  (x/c>7)is  not  large,  and  the  lift  enhancement  is  mainly  produced 
by  the  overall  aerodynamic  change  on  the  main  airfoil  induced  by  SETE.  Comparisons  of  the  calculated  and 
measured  lift  coefficients  are  shown  in  Fig.  18  as  a  function  of  AoA.  The  calculated  results  for  the  2D  airfoil  are 
converted  to  the  3D  wing  results  using  Eq.  (7).  The  calculated  and  experimental  results  are  consistent  until  the 
calculated  lift  coefficients  level  at  higher  AoAs.  As  indicated  in  Fig.  13,  the  calculated  increment  in  the  lift 
coefficient  exhibits  the  same  trend  as  the  experimental  and  theoretical  results  as  a  function  of  the  deflection  angle  of 
SETE.  Figure  19  shows  velocity  magnitude  contours  around  the  NACA0012  airfoil  with  zero  AoA  for  0°  and  5° 
deflection  angles.  The  asymmetry  of  the  flow  field  is  induced  by  SETE  and  the  wake  is  turned  downward, 
indicating  a  deflected  momentum  stream  tube  and  the  generation  of  additional  lift.  The  wake  structure  is  not 
appreciably  altered  although  it  is  deflected;  this  implies  that  the  parasite  drag  is  not  significantly  affected. 

Further,  a  10%c  SETE  with  a  5°  deflection  angle  is  compared  with  a  conventional  plain  flap  that  is  deflected  at  the 
90%  chord  of  a  NACA0012  airfoil.  Figure  20  shows  the  computed  lift  and  drag  coefficients  as  a  function  of  the 
deflection  angle  for  the  SETE  and  corresponding  plain  flap.  Compared  with  the  plain  flap,  the  lift  enhancement  by 
the  SETE  is  larger  while  the  drag  coefficient  is  smaller.  Here  the  lift  and  drag  coefficient  for  the  SETE  is  defined 
based  on  the  main  airfoil  chord  c.  Even  when  these  coefficients  are  based  on  1.1c  for  the  SETE,  the  above  results 
are  still  qualitatively  valid.  The  pressure  coefficient  distributions  on  the  main  airfoil  for  the  SETE  (5°  deflection), 
plain  flap  (5°  deflection)  and  Gurney  flap  are  shown  in  Fig.  21,  indicating  that  the  larger  pressure  difference  on  the 

main  airfoil  is  induced  by  the  SETE.  The  Gurney  flap  height  is  0.1csin(  5""  )  »  0.87%c  ,  which  is  the  same  as  the 
maximum  vertical  displacement  of  the  SETE.  Figures  22  and  23  show  the  lift  coefficient  as  a  function  of  AoA  and 
the  drag  polar  for  the  SETE,  plain  flap  and  Gurney  flap.  Among  these  trailing  edge  devices,  the  SETE  generates  the 
largest  lift  increase  at  a  smallest  drag  penalty.  It  should  be  emphasized  here  that  the  absolute  accuracy  of  the  drag 
predicted  by  CFD  is  not  critically  examined  although  the  computed  drag  coefficient  for  the  NACA0012  airfoil  is 
close  to  Cjy  =0.02  -  0.03  given  by  measurements  for  the  corresponding  airfoil  section.  Therefore,  the  above 
comparison  in  the  drag  polar  is  valid  only  in  a  relative  sense.  Furthermore,  for  the  2D  airfoil,  the  drag  polar  is  not 
parabolic  since  no  induced  drag  is  included. 
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Figure  17.  Calculated  pressure  coefficient 
distributions  on  the  NACA0012  airfoil  for  the  SETE 
deflection  angles  of  0  and  5  degrees  at  zero  AoA. 


Figure  18.  Comparisons  of  the  calculated  and 
measured  lift  coefficients  as  a  function  of  AoA  for  the 
deflection  angles  of  0  and  5  degrees. 
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(a)  (b) 

Figure  19.  Velocity  magnitude  fields  of  the  NACA0012  airfoil  at  zero  AoA,  (a)  zero  deflection  of  SETE,  (b)  5 
degrees  deflection  of  SETE. 


Figure  20.  Comparison  of  the  computed  lift  and  drag  coefficients  between  SETE  and  conventional  flap  for  a 
NACA0012  airfoil  at  zero  AoA. 


Figure  21.  Comparison  of  the  computed  pressure  coefficient  between  SETE,  conventional  flap  and  Gurney 
flap  for  a  NACA0012  airfoil  at  zero  AoA. 
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Figure  22.  The  computed  lift  coefficient  as  a  function 
AoA  for  SETE  (5-deg  deflection),  Gurney  flap  and 
conventional  flap  (5-deg  deflection)  for  a  NACA0012 
airfoil. 

Conclusions 
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Figure  23.  The  computed  drag  polar  for  SETE  (5-deg 
deflection),  Gurney  flap  and  conventional  flap  (5-deg 
deflection)  for  a  NACA0012  airfoil. 


Experiments  and  calculations  presented  show  that  a  static  extended  trailing  edge  (SETE)  attached  to  a  NACA0012 
airfoil  model  is  able  to  enhance  the  lift  while  the  zero-lift  drag  is  not  significantly  increased.  The  mechanism  of  lift 
enhancement  by  SETE  is  the  camber  effect,  which  is  the  same  as  that  for  other  high-lift  devices  like  Gurney  flap  and 
conventional  flap.  However,  compared  with  Gurney  flap  and  conventional  flap,  SETE  generates  a  larger  lift 
increase  at  a  smaller  drag  penalty  since  it  is  embedded  in  the  wake  of  the  main  airfoil.  Therefore,  SETE  has  a 
promising  potential  for  improving  the  efficiency  of  cruise  flight.  The  benefit  margin  of  SETE  for  cruise  flight  is 
evaluated  in  comparison  with  Gurney  flap,  and  the  feasibility  of  SETE  for  lift  enhancement  in  cruise  flight  is 
demonstrated  at  small  angles  of  attack  and  deflection  angels  of  SETE.  The  mechanical  simplicity  of  SETE  allows 
direct  application  to  aircraft  without  changing  the  basic  aerodynamic  configuration.  Eurthermore,  the  deformation 
and  pitching  moment  of  SETE  due  to  aerodynamic  loading  are  estimated,  and  actuating  SETE  is  feasible  since  the 
pitching  moment  on  it  is  small.  Euture  research  will  focus  on  SETE  actuation  using  smart  material  like 
piezoactuators  and  shape  memory  alloys  for  steady  and  unsteady  flow  control. 
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Flexible  Extended  Trailing  Edge 

—  A  Biologically-Inspired  Concept 


Objectives 


Main  Airfoil 


Extended  Trailing  Edge 


•  Understand  the  aerodynamics  characteristics  of  static 
extended  traiiing  edge  (SETE)  via  integrating  CFD  and 
experiments  (EXP). 

•  Explore  the  feasibility  of  improving  the  aerodynamic 
performance  of  airfoils  and  wings  using  flexible 
extended  trailing  edge 


Steady  and  Unsteady  Aerodynamic  Aspects 

SETE  for  Lift  Enhancement  and  Drag  Reduction 
in  Cruise  Flight  (small  Ao A): 


Main  Airfoil 


Extended  Trailing  Edge 


Dynamic  flexible  extended  trailing  edge  for 
Separation  Control  (high  AoA): 


feedback 


flexible  TE 


Technical  Approaches  — 

Combination  of  Experimental,  Computational  and  Theoretical  Methods 

•  Experiments: 

Quantities:  Integrated  forces,  pressure,  velocity  fields, 

shear  stress,  shape,  kinematics,  strain 
Integrated  Techniques:  Balance,  pressure  transducers, 

PIV,  high-speed  stereo  videogrammetry 
shear-sensitive  LC,  TSP,  PSP,  strain  gauges, 
oil  film  skin  friction  meter,  etc. 

Responsive  Skin:  Sensing  and  actuating  polymer 

•  CFD:  PANS  and  LES  for  flows  coupled  with  finite  element 

code  for  flexible  structure 

•  Theoretical  Study:  Adapted  thin-airfoil  theory,  unsteady 

extension  coupled  with  thin-plate  dynamics, 
interaction  between  BL  and  wake  via  global  stability 
analysis 


Teaming  — 

Combination  of  Experimental,  Computational  and  Theoretical  Methods 
Organization:  ^ 

Western  Michigan  Univ.  -  CFD  &  EXP  *. 

NASA  Langiey  Research  Center  -  EXP 

Personnel:  **  ^ . 

Faculty  (2),  Research  Scientist  (1),  Post-Doc  (1),  GS-Doc  (1), 
GS-MS(1),  UG  students(2). 

Faculty  One 


Post  Doc 

^SS-MS 
UG  Students 
Research  Scientist 


Communications: 


Meetings,  Emails,  Calls,  Site  Visits,  Hall-Way  Conversations. 
Immersive  Integration. 


Computational  Fluid  Dynamics  Calculations 

Solver: 

Commercial  Code:  RANS  using  FLUENT,  CFX 
In-House  Code:  3D  RANS  solver  with  immerse 
boundary  method,  Incompressible,  Second-Order, 
Turbulence  Modeling. 

Meshing: 

Structured  and  unstructured  using  ICEM  and  in-house 
Grid  Independence.  Grid  Convergence. 


Validation: 

NACA0012  Laminar  and  Turbuient. 


Aerodynamic  Force  Measurements 

in  Wind  Tunnels 

NACA0012  Model  with  SETE  (provided  by  NASA  LaRC) 


Advanced  Design  Wind  Tunnel: 

Speed:  6-73  m/s 
Test  section:  4  by  3  feet 
Tu:  0.1 -0.4% 

Six-component  balance 
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L/D  of  NACA0012  with  SETE  EXP 


•  Zero-lift  drag  and  Oswald  efficiency  remain  unchanged 
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Comparison  between  SETE  and  Conventional  Flap 
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Comparison  of  SETE,  Gurney,  and  Conventional  Flap 


Main  Airfoil 


Extended  Trailing  Edge 


Main  Airfoil 


Gurney  Flap 


CFD 


Main  Airfoil 


Flap 


Comparison  of  SETE,  Gurney,  and  Conventional  Flap 

CFD 


0  0.2  0.4  0.6  0.8  1 

x/c 


CL  Increment 


Thin- Airfoil  Theoretical  Interpretation  for  SETE 
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Concluding  Remarks 

Static  extended  trailing  edge  (SETE)  attached  to  a 
NACA0012  airfoil  model  is  able  to  enhance  the  lift 
while  the  zero-lift  drag  is  not  significantly  increased. 

Camber  effect  caused  by  the  SETE  on  the  main  airfoil. 

Compared  with  Gurney  and  conventional  flaps,  SETE 
generates  larger  lift  with  smaller  drag  penalty-  ideal  for 
cruise  flights. 

SETE  is  mechanically  simple  to  implement  to  existing 
platforms. 


Responsive  Skins  for 
Flexible  Extended  Trailing  Edge  Control 


•  Birds  utilize  their  wing  flexibility  for  their 
different  flight  regimes. 

•  Airfoil  section  with  an  extended  trailing 
edge  embedded  with  MEMS  microphones. 

•  The  responsive  skin  (trailing  edge  embedded 
with  MEMS  microphones)  will  sense 

and  react  to  changes  in  surroundings. 


Pressure  sensors 


MEMS  Microphones 
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Flexible  MEMS  Array 

(Pressure,  Humidity,  Temperature,  Mierophone,  2-axis  Aeeelerometer) 


Two-axis  accelerometei 


MEMS  microphones 


SOI  (Silicon  on  Insulator) 
High  Temp.  Press.  &  Tem 
Up  to  240  C 


Absolute  Pressure 

Pressure, 

Humidity,  Temp. 


MEMS  sensors  on  flexible  skin  have  been  designed 
and  will  be  characterized. 
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Concluding  Remarks 


Integrating  CFD  and  experiments  capabilities  in  the 
research  process  provide  countless  benefits: 

-  Scholarship 

-  Cost 

-  Time  to  Market 

-  Team  Leadership 

-  Education 
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